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Isokinetic Sampling 

BACKGROUND OF THE DISCLOSURE 

[0001] The present invention relates to isokinetic 
sampling. It is particularly, but not exclusively, 
related to methods and apparatuses for carrying out 
isokinetic sampling, and in particular isokinetic 
sampling in oilfield applications. 

[0002] In oilfield applications, as in many other 
fields, it is important to be able to analyse the 
composition and properties of a multi-phase fluid stream, 
such as a gas-condensate stream or a wet-gas stream. In 
particular it is desired to know the guantity of the 
various phases of the flow, for example by knowing the 
gas-oil ratio (GOR) or the condensate-gas ratio (CGR) . 
It is also desired to know the properties of the various 
phases, for example their pressure, volume and 
temperature (PVT) relationships. 

[0003] For example, the prediction of fluid phase 
behavior and reservoir simulation models based on 
eguations-of-state (EOS) typically reguires high guality 
PVT data. 

[0004] PVT data are also often needed for flow 
assurance in wells and transport lines. 

[0005] Furthermore, PVT data often have a significant 
impact on processing facility designs and specifications, 
and therefore on the profitability of gas-oil field. 
[0006] However, to accurately determine the PVT 
properties and composition of reservoir fluid, 
representative fluid samples are reguired. Common 
analysis sampling procedures suffer from deficiencies in 
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either or both of their accuracy and their ability to 
cope with high flow rates. 

[0007] For example, a high-rate gas condensate well 
producing 5-100 MMscf (Millions standard cubic feet) /day 
with a gas-oil ratio of 3k-100k scf/bbl (CGR from 0 to 
200 bbl/MMscf) (bbl = barrel of liguid; 1 m 3 = 6.29 bbl), 
can often overwhelm a test-separator, causing liguid 
carry-over in the separator gas outlet line, thereby 
providing poor measurement of the GOR and non- 
representative PVT samples (if the carry-over is not 
measured and corrected for) . This normally also results 
in poor recombination ratios. Gas condensate wells are 
particularly problematic for sampling as large volumes 
of gas are associated with only small volumes of liguid, 
and the phase behavior will be highly sensitive to the 
guantity and composition of the liguid phase. 
[0008] Wellhead sampling is therefore considered to be 
the only practical method of obtaining reliable test data 
for such applications. Isokinetic sampling at a wellhead 
of a multiphase fluid is desirable since, if achieved, it 
means that the sample of the fluid is at the same 
pressure, temperature and velocity as the main flow 
stream, and therefore will have identical properties to 
that main flow stream. 

[0009] A known wellhead sampling system is shown 
schematically in Figure 1, which is taken from Dybdahl 
and Hjermstad "A systematic approach to sampling during 
well testing", Society of Petroleum Engineers (SPE) paper 
69427, SPE Latin America & Caribbean Petroleum 
Engineering Conference, Buenos, Argentina, 25-28 March 
2001 (see also the paper by Henk et al . "Testing of gas 
condensate reservoirs - sampling, test design and 
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analysis", SPE paper 68668, SPE Asia Pacific Oil & Gas 
Conference, Jakarta, Indonesia, 17-19 April 2001) . Such 
a system is, for example, produced by Petrotech under 
their IsoSplit (TM) wellhead sampling and separator 
sampling services. 

[0010] The wellhead sampling system of Figure 1 
consists of two parts. 

[0011] The first part is a wellhead mixing/sampling 
manifold 1010, which is positioned upstream of the choke 
manifold 1030, and has a f ocusing-action static mixing 
device 1015 (see GB 2301297 & US 5894080; a similar 
mixing device is shown in Figure 2a) for collecting 
liquid at the wall of the pipe and distributing it into 
the main flow stream. A sampling probe 1020 is radially 
inserted downstream of the mixer 1015 for withdrawing a 
sample of the main flow stream. 

[0012] The second part is an on-site laboratory 1050 
(also called a "Mini-lab"), which processes the sampled 
stream to measure the condensate-gas ratio (CGR) . PVT 
samples for compositional characterization are taken. 
[0013] As described in SPE paper 68668, the liquid 
distribution across the flow pipe of Figure 1 will vary, 
and either a mixing device or a traversing probe has to 
be used in order to have a split stream sample with the 
correct gas condensate ratio. As described in SPE paper 
69427, the probe traversing enables sampling at multiple 
points across the pipe and the subsequent correction for 
the droplet and velocity distribution. It is stated in 
SPE 68668 that the liquid collected from the wall by the 
f ocusing-action mixer 1015 is distributed homogeneously 
into the flowing gas stream. However, our own gas-liquid 
ratio (GLR) profile measurement shown in Figure 2b, for a 
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mixer design (shown in Figure 2a) similar to the 
focussing-action mixer 1015, indicates a liguid-rich jet 
at the central region of the sampling cross section. The 
GLR measured at such points will be a gross under- 
estimate of the pipe average value. Figure 2b shows a - 
70% error in the GLR at the pipe centre, with respect to 
the pipe reference GLR. 

[0014] SPE paper 68668 also states that sample 
withdrawal through the sampling probe 20 from the main 
stream is performed at an isokinetic sampling rate, i.e. 
at the same linear velocity as the wellhead fluid stream, 
which is necessary to sample the gas and liguid phase in 
the correct ratio. This isokinetic sampling rate is said 
to be calculated from the probe size and the gas flow 
rate measured at the test separator. 
[0015] This approach has at least the following 
disadvantages : 

1) A test separator is needed for wellhead 
sampling operations, which is bulky and 
expensive (Figure 1). 

2) The test separator can be overwhelmed by 
high-rate gas wells, leading to liguid 
carry-over in the separator gas outlet line. 
This can result in an inaccurate main-stream 
gas rate measurement used to set the 
isokinetic sampling rate. 

3) The main-flow gas rate can fluctuate, for 
example due to changes in the liguid 
loading. This unsteadiness in the main flow 
rate can cause the separator gas rate 
measurement to be unreliable (see SPE paper 
76 76 6, by B.C. Theuveny and P. Mehdizadeh 
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"Multiphase flowmeter application for well 



and fiscal allocation", Alaska, 20-22 May 



2002), leading to unreliable isokinetic 



sampling. 



4) 



Isokinetic sampling rate control based on 



measured main-stream gas flow rates can lead 
to errors when liguid loading becomes 
relatively high. 



[0016] WO00/49370 describes a method for measuring the 
liguid and gas flow rates of a multi-phase fluid stream 
including withdrawing a portion of the overall flow rate 
under "isokinetic conditions". The sampled portion forms 
approximately 5% to 15% of the total flow rate. 
[0017] The methods and apparatus disclosed in 
WO00/49370 do not independently measure the main flow 
stream, but instead derive measurements of the main flow 
stream from the effect of the sampling operation on the 
flow stream. Furthermore, although the methods described 
are claimed to be "isokinetic", this term is applied to 
all sampling in which the ratio of the sampled flow rate 
to the calculated total flow rate is within 20% of the 
ratio of the cross sectional area of the sampling probe 
to the cross sectional area of the flow stream. 
Adjustment of the sampling conditions is only made if the 
ratio falls outside this range. 

[0018] The method proposed in WO00/49370 for measuring 
the total flow rate (CM of the fluid stream is given in 
eguation (5) of that application, with Q 1 being calculated 



as the sole unknown in the eguation — --\ — , where q 

Ap 2 {Q l -q) 



is the measured sampled flow rate and Ap 2 and Ap 2 are 
measured pressure differences of the main flow stream 
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with and without sampling operating. This assumes that 
the density of the sampled flow q is the same as the 
density of the total flow Q t . 

_ ^A Pl IAp 2 
[0019] Thus Q x is calculated as U\ - Q r- — ^ . 

yjAp 1 / Ap 2 -1 

From this equation, the error 5Q 2 in Qj can be determined 
« lfj{fi + ifj=lf- " here P = V^- For P ,> 1, 
= — . However, as J Ap, I Ap, — > 1 , — >°°. Therefore, in 

order to avoid excessive errors in the measurement of the 
total flow rate Q 1 when using the method of WO00/49370 it 
is important that the difference between Ap 2 and Ap 2 is 
large, and consequently that the flow rate of the sampled 
flow is a substantial proportion (i.e. more than 5%) of 
the total flow rate. When such sampling proportions are 
used in high flow situations, for example in a high-rate 
gas condensate wells, they can overwhelm the separator, 
or alternatively require a larger, and therefore more 
expensive separator. 

BRIEF SUMMARY OF THE DISCLOSURE 
[0020] The present invention is at least partly based 
on the realisation that isokinetic sampling from a fluid 
stream can be achieved by monitoring the ratio of the 
flow rate of the sampled flow to the total flow rate of 
the main fluid stream and adjusting the sampling rate 
accordingly. Contrast this with WO00/49370 which 
provides no suggestion that any repeated or continuous 
adjustment to achieve isokinetic sampling should be 
carried out, and indeed assume that sampling is carried 
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out under isokinetic conditions in order to measure the 
total flow rate. 

[0021] Thus, a first aspect of the present invention 
provides a method of sampling a multi-phase fluid stream 
comprising the steps of: 

sampling, with a sampling probe, a portion of the 
fluid stream; 

measuring the flow rate Q s of said sampled portion; 

and 

measuring, independently of the sampling step, the 
total flow rate Q of the fluid stream, 

wherein the flow rate of the sampled portion is 
controlled according to the ratio of the flow rate of the 
sampled portion to the flow rate of the fluid stream, in 
order to obtain substantially isokinetic sampling of the 
fluid stream. 

[0022] Preferably the flow rate of the sampled portion 

Q A b 

is controlled such that the eguation — = pmbe is 

Q \, PC 

satisfied to within 10%, more preferably to within 5%, 
where A probe and A plpe are the respective cross-sectional 
areas of the sampling probe opening and the fluid stream 
at the sampling point. 

[0023] The flow rates Q and Q s may be mass flow rates, 
and the flow rate Q s of the sampled portion of the fluid 
stream is preferably less than 0.2% of the total flow 
rate Q, more preferably less than 0.05% of the total flow 
rate Q. 

[0024] The method may be used in high-rate gas 
condensate wells, where the gas flow rate is at least 5 
MMscf/day. In wells with higher flow rates, the flow 
rate of the sampled portion of the fluid stream is 
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particularly preferred to be a small proportion of the 
total flow rate so that downstream analytical processes 
on the sampled portion are not overwhelmed. An advantage 
of the present invention is that it can be used to sample 
relatively small proportions of the total flow. 
[0025] Preferably the method further includes the step 
of separating said sampled portion into a sampled liguid 
flow and a sampled gas flow, and said step of measuring 
the sampled flow rate includes measuring the flow rate of 
the sampled liguid flow and the flow rate of the sampled 
gas flow. Either or both of these sampling flows may be 
measured using a Coriolis flow meter. This allows more 
accurate measurement of the flow rates of the different 
phases of the fluid stream. A Coriolis flow meter also 
allows the density of the flow to be measured. 
[0026] In certain embodiments the step of measuring 
the total flow rate includes measuring the density of the 
fluid stream. The step of measuring the total flow rate 
may include measuring a differential pressure, which may, 
for example, be measured using a Venturi meter, or be an 
impact pressure measured using a Pitot tube. 
[0027] Preferably the step of measuring the total flow 
rate is carried out continuously or repeatedly during the 
sampling process. 

[0028] Preferably the step of controlling the flow 
rate of the sampled portion is carried out continuously 
or repeatedly during the sampling process. 
[0029] The continuous or repetitive measuring of the 
total flow rate and controlling of the flow rate of the 
sampled portion allows correction of the flow rate of 
sampled portion in response to changes in the mass flow 
rate of the main flow stream. 
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[0030] In this respect, the term "continuous" also 
includes repeated discrete applications of the measuring 
and controlling steps above at regular intervals, these 
intervals being shorter than the time scale over which 
significant changes in the total mass flow rate of the 
main stream occur. 

[0031] The above method is particularly applicable 
where the fluid stream is a feed from a well-head, and 
may be carried out at a well-head site. 

[0032] Preferably the method further includes, prior 
to sampling the fluid stream, the step of conditioning 
the fluid stream. 

[0033] Another aspect of the present invention 
provides a sampling apparatus for a multiphase fluid 
stream which can carry out the method of the above 
aspect . 

[0034] According to this aspect, the present invention 
provides a sampling apparatus for sampling a multi-phase 
fluid stream, comprising: 

means for measuring the total flow rate Q of the 
fluid stream; 

a sampling probe for sampling a portion of said 
fluid stream; 

a meter for measuring the flow rate Q s of said 
sampled portion; and 

a controller adapted to control the flow rate of the 
sampled portion according to the ratio of the flow rate 
of the sampled portion to the total flow rate of the 
fluid stream, in order to obtain substantially isokinetic 
sampling of the fluid stream. 

[0035] Preferably the controller is adapted to control 
the flow rate of the sampled portion such that the 



57.0503 U.S. PCT 



10 

equation — - = proh " is satisfied to within 10%, more 

Q \ ve 

preferably to within 5%, where A probe and A pipe are the 
respective cross-sectional areas of the sampling probe 
opening and the fluid stream at the sampling point. 
[0036] Flow rates Q and Q s may be mass flow rates, in 
which case the meter is a mass flow rate meter. 
[0037] Preferably the cross-sectional area of the 
sampling probe opening is less than 0.2% of the cross- 
sectional area of the fluid stream at the sampling point, 
and thus the flow rate of the sampled portion is a 
similar fraction of the total flow rate of the fluid 
stream. More preferably the cross-sectional area of the 
sampling probe opening is less than 0.05% of the cross- 
sectional area of the fluid stream at the sampling point. 
[0038] The sampling apparatus may form part of a high- 
rate gas condensate well wherein the gas flow rate is at 
least 5 MMscf/day. 

[0039] Preferably the apparatus further includes a 
separator for separating the sampled portion into a 
sampled liquid flow and a sampled gas flow; a meter for 
measuring the flow rate of the sampled liquid flow; and a 
meter for measuring the flow rate of the sampled gas 
flow. Either or both of these meters may be Coriolis 
flow meters . 

[0040] In some embodiments the means for measuring the 
total flow rate may include a densitometer for measuring 
the density of the fluid stream. The means for measuring 
the total flow rate may additionally or alternatively 
include a differential pressure sensor, such as that for 
use with a Venturi meter or a Pitot tube. 
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[0041] Preferably the means for measuring the total 
flow rate operates continuously or repeatedly during the 
sampling process. 

[0042] Preferably the controller controls the flow 
rate of the sampled portion continuously or repeatedly 
during the sampling process. 

[0043] Preferably the sampling probe is made from 
damage resistant material such as tungsten carbide or 
cobalt alloy, so that the cross sectional area of the 
sampling probe opening does not vary during the sampling 
process . 

[0044] The sampling probe is preferably located in the 
centre of the pipe section. 

[0045] Preferably the sampling apparatus further 
includes a flow conditioner located upstream of the 
sampling probe . 

[0046] The fluid stream may be a feed from a well- 
head . 

[0047] In order to achieve representative isokinetic 
sampling of a multi-phase stream using sampling at a 
single sampling point, it is preferable that, at the 
designated sampling cross-section of the pipe, the cross- 
pipe GLR profile is uniform. 

[0048] Previously flow-focussing mixers have been used 
to try to provide such a uniform GLR profile. Examples 
of such mixers are plug-valve mixers, for example as 
shown in GB 2301297 and US 5894080. A similar mixer to 
that in GB 2301297 is shown schematically in Figure 2a. 
Where similar reference numerals are used in Figure 2a, 
they refer to features already described in association 
with Figure 1 . 
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[0049] However, use of flow-focussing mixers, such as 
the mixer 1015 shown in Figure 2b, can result in too much 
liguid being entrained into the central pipe region as a 
result of the focussing action of the mixers. At low 
pressures there may be excessive liguid jetting, leading 
to very non-uniform cross-pipe GLR and a large GLR error. 
Figure 2b shows the GLR cross-pipe profile for air/water 
flows with an eguivalent CGR of -60 bbl/MMscf gas (gas 
rate g g ~ 1400 m 3 /h or about 1.2 MMscf / day) , and shows 
that the measured GLR has an error of approximately -70% 
from the reference value at the pipe centre. 
[0050] It is anticipated that as the wellhead pressure 
increases, the GLR cross-pipe non-uniformity will improve 
as a results of the reduction in the gas-liguid density 
contrast, and hence in the gas-liquid velocity slip. 
However, the adverse effect of liquid jetting will still 
remain, even at high wellhead pressures. Similar results 
are expected to be seen with the plug-valve type mixers 
referred to above. 

[0051] Therefore a further aspect of the present 
invention provides a flow conditioner for a wellhead 
sampling system, the flow conditioner being a pipe 
section including: a swirl inducing section; a flow 
straightener ; and an orifice plate. 

[0052] To achieve accurate isokinetic sampling at 
wellhead at a single sampling point, it is preferable to 
have a uniform cross-pipe GLR profile at the designated 
sampling cross-section. For annular or annular-mist 
flows prevailing in a wellhead sampling system (i.e. not 
fully developed vertical upward or downward flow) , this 
can be achieved by removing liguid film from the pipe 
wall and by entraining it homogeneously in the gas core. 
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It is also preferable that liquid drop re-deposition onto 
the pipe wall should be negligible at the sampling cross- 
section. Drop size is preferably small thereby 
minimizing gas-drop velocity slip and the rate of drop 
deposition. These factors will lead to a more 
representative sampling of gas and liquid flow rates, and 
hence a better measure of the GLR. 

[0053] The creation of the initial swirl, upstream of 
the orifice plate, preferably creates a uniformly 
circumf erentially distributed liquid film on the wall of 
the pipe section. The orifice plate sheds any such 
liquid film on the pipe wall back into the gas flow. 
After the initial formation of a liquid rich core 
downstream of the orifice plate, turbulence causes good 
mixing of the liquid and gas phases. Straightening the 
flow of the gas core prior to the orifice plate reduces 
the amount of liquid that is redeposited onto the pipe 
wall after shedding by the orifice plate. Therefore a 
relatively uniform cross-pipe GLR can be achieved, making 
it less important where on the pipe diameter in the 
sampling cross section a sampling probe is placed. 
[0054] The swirl inducing section may include a 
tangential input to the pipe section. Alternatively the 
swirl inducing section may include a helical insert 
inside the pipe section. 

[0055] The flow straightener is preferably located in 
the centre of the pipe section and such that there is an 
annular gap between the flow straightener and the wall of 
the pipe section, thereby straightening the gas core of 
the fluid stream whilst not affecting any liquid film 
located on the wall of the pipe. 
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[0056] The orifice plate may have a diameter of its 
open section which is between 0.85 and 0.95 of the local 
diameter of the pipe section, preferably between 0.9 and 
0.95 of the local diameter of the pipe section (also 
referred to as the diameter ratio f3) . 

[0057] The pipe section may narrowed at the orifice 
plate. The narrowed pipe section gives a higher gas flow 
momentum and higher gas Weber number. 

[0058] To improve CGR measurement accuracy of lean gas 
condensate wells (CGR = 0 to 50 bbl/MMscf), a dual- 
orifice mixer arrangement is preferably used for reducing 
liguid rate under-sampling at low liguid loadings. 
[0059] In one example of such an arrangement, the 
first (i.e. upstream) orifice plate has a diameter ratio 
such as ($1=0.9, and the second orifice plate has a 
diameter ratio such as /3 2 =0.95. The two orifice plates 
may be axially separated by ID to 3D, where D is the 
basic diameter of the pipe at the orifice plates. 
[0060] The first, upstream orifice plate (diameter 
ratio /Si) is intended to shed as much as liquid to the 
pipe centre as possible (with a relatively small /3 - 
ratio) ; whilst the second orifice plate is used to shed 
any relatively thin liquid film that might be caused by 
the re-deposition of the liquid drops shed by the first 
orifice plate, after travelling a few pipe diameters. 
[0061] A sampling probe may be used to sample the 
fluid stream, and is preferably located downstream of the 
orifice plate (s). More preferably the sampling probe is 
located at a distance between ID and 3D downstream of the 
downstream orifice plate, where D is the diameter of the 
pipe section at the orifice plate (not the diameter of 
the orifice plate itself) . 
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[0062] Preferably the sampling probe is located in the 
centre of the pipe section. 

[0063] Another aspect of the present invention 
provides a method for conditioning the flow of a multi- 
phase stream comprising the steps of successively: 

forming a liguid film at the wall of a pipe section 
which carries the fluid stream; 

straightening the flow of at least a core of the 
fluid stream; and 

shedding the liguid film into said core. 

Preferably the liguid film is substantially 
uniformly circumf erentially distributed on the pipe wall. 
[0064] The step of shedding the liguid film may use at 
least one orifice plate. The or each orifice plate may 
have an open diameter of 0.85 to 0.95 of the diameter of 
the pipe section in which it is placed, preferably 
between 0.9 and 0.95 of the diameter of the pipe section. 
The pipe section may narrow at at least the orifice 
plate . 

[0065] The step of forming the liguid film may include 
causing the fluid stream to swirl about the axis of the 
pipe section. Such rotary motion may be achieved by, for 
example, tangentially introducing the fluid stream to the 
pipe section, or by passing the fluid stream through a 
helical insert inside the pipe section. 

[0066] Preferably the step of straightening the core 
of the fluid stream does not straighten a portion of the 
fluid stream adjacent to the walls of the pipe, and 
therefore does not affect the liguid film. 
[0067] The fluid stream may be sampled using a 
sampling probe, the sampling preferably taking place 
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downstream of the orifice plate (s). More preferably the 
sampling probe is located at a distance between ID and 3D 
downstream of the downstream orifice plate, where D is 
the diameter of the pipe section at the orifice plate 
(not the diameter of the orifice plate itself) . 
[0068] Preferably the sampling probe is located in the 
centre of the pipe section. 

[0069] The flow conditioner or method of the above two 
aspects may be used with the sampling methods and 
apparatuses of the earlier aspects to improve the 
accuracy of the isokinetic sampling. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0070] Embodiments of the present invention will now 
be described in relation to the accompanying drawings, in 
which : 

[0071] Figure 1 is a schematic illustration of a known 
wellhead sampling system, and has already been described; 
[0072] Figure 2a is a schematic illustration of a 
flow-focusing mixer similar to that shown in GB 2301297 
the prior art and has already been described; 
[0073] Figure 2b shows the GLR cross-pipe profiles 
downstream of the mixer of Figure 2a and downstream of a 
flow conditioner according to an embodiment of the 
present invention; 

[0074] Figure 3 is a flow chart showing a method 
according to an embodiment of the present invention; 
[0075] Figure 4 is a schematic illustration of an 
isokinetic sampling system according to an embodiment of 
the present invention; 

[0076] Figure 5 is a schematic illustration of an 
isokinetic sampling system according to another 
embodiment of the present invention; 
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[0077] Figure 6 is a schematic illustration of an 
isokinetic sampling system according to a further 
embodiment of the present invention; 

[0078] Figure 7 is a schematic diagram of a flow 
conditioner according to an embodiment of the invention; 
[0079] Figure 8 is a schematic diagram of a flow 
conditioner according to an alternative embodiment of the 
invention; and 

[0080] Figure 9 is a schematic diagram of a flow 
straightener which may be used in the flow conditioners 
shown in Figures 7 and 8. 

DETAILED DESCRIPTION 

[0081] Elements of a method of sampling according to 
an embodiment of the present invention are shown in the 
flow chart of Figure 3. 

[0082] A portion of a fluid stream is sampled, and the 
flow rate Q s of that sampled portion measured. The total 
flow rate Q of the fluid stream is also independently 
measured, and the ratio of the two measured flow rates 
calculated. Based on this ratio, the sampling flow rate 
(i.e. the proportion of the fluid stream which is 
sampled) is adjusted in order to substantially obtain 
isokinetic sampling conditions. This process is 
preferably either repeated or carried out continuously 
throughout the sampling process. 
[0083] The sampled portion may be analysed 
concurrently with the flow rate measurements and 
adjustments, for example to provide the GOR, the WLR and 
samples for PVT analysis. This analysis may also use the 
flow rate measurement of the sampled portion. 
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[0084] Figure 4 shows an isokinetic sampling system 
according to an embodiment of the present invention. 
[0085] An isokinetic sampling model is used to 
calculate the target sampling flow rate Q so , according to 
the main stream flow rate Q e ff- One example of such an 

Qso A pro be 

isokinetic sampling model is = for substantially 

Qeff A pipe 

uniform mixed flow. 

[0086] In contrast to the method described in relation 
to Figure 1 above, in the present embodiment, the main 
stream flow rate Q eff is measured on-line in the pipe by a 
main pipe flow meter 40. This may be, for example, a 
Venturi or an orifice meter. The measured Q eff is then 
fed into the sampling model so that Q so is updated as Q e ff 
varies . 

[0087] When the on/off valve 35, allowing sampling to 
take place, the actual sampling flow rate Q s is measured 
by a sampling flow meter 30 positioned in the sampling 
line. The actual sampling flow rate Q s and the target 
sampling flow rate Q so are compared in comparator 33, 
which creates a flow control signal 31, on the basis of 
which control valve 32 is adjusted by the actuator 34 so 
as to keep Q s = Q so . 

[0088] Alternatively flow control signal 31 may be 
created as a result of a direct comparison between the 
main stream flow rate Q eff measured by the main pipe flow 
meter 40 and the sampling flow rate Q s measured by the 
sampling flow meter . 

[0089] The main pipe flow meter 40 of the embodiment 
of Figure 4 may be positioned upstream or downstream of a 
mixer or other flow conditioner (not shown) . In certain 
embodiments, the main pipe flow meter will also serve as 
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a good gas-liquid mixer, and so a separate flow 
conditioner may not be required. 

[0090] However, use of a flow conditioner upstream of 
a flow meter 40 can serve to eliminate or substantially 
reduce the gas-liquid velocity slippage for the flow 
passing through the main pipe flow meter 40, thereby 
allowing a homogenous flow model (e.g. standard single- 
phase Venturi equations) rather than a slip-model (e.g. a 
2-layer, separated model for annular or annular-mist 
flows) to be used to calculate the main stream flow rate. 
The homogenous mixing state is assumed to be maintained 
over some distance downstream of the main pipe flow meter 
40, enabling the sampling probe to perform the desired 
isokinetic sampling. 

[0091] Preferred embodiments of separate mixers or 
flow conditioners are described below. 

[0092] Another embodiment of the present invention is 
shown schematically in Figure 5. An embodiment of the 
method of the present invention is also described in 
relation to this Figure. Again components which have 
already been discussed in relation to Figures 1 and 4 
have the same reference numbers. 

[0093] This embodiment has a flow meter 40 which 
combines a Venturi pressure drop measurement 41 with a 
gamma-ray densitometer 42. An example of such a meter is 
the PhaseTester/PhaseWatcher (TM) Vx Venturi meter 
described in SPE paper 77405 (by S. Jayawardane and B.C. 
Theuveny "PVT Sampling With Multiphase Flowmeters- 
Theoretical Justifications and Field Limitations", SPE 
Annual Technical Conference & Exhibition, San Antonio, 
Texas, 29 September to 2 October 2002) . The pressure and 
temperature (p, T) of main flow stream 11 can also be 
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measured at this point. The Vx flow meter has been shown 
(see SPE paper 76766) to give better multiphase flow rate 
measurement than a test separator, especially at high gas 
rate and/or unsteady flow conditions. 

[0094] This embodiment uses a total mass flow rate 
measurement Q = K^Ap.p H , obtained from the combined 
measurements of high-energy mixture density (p Hr 
calculated to about 10 kg/m 3 accuracy by the Vx nuclear 
source) and Ap from the Venturi measurement (calculated 
to within 1% accuracy) to set the extraction rate for 
isokinetic sampling, where K is a calibration factor 
depending on Venturi throat diameter, inlet diameter and 
the discharge coefficient; see, for example, eq.5 in SPE 

AO i[a(Ap) A/? h 1 

paper 77405. The mass flow rate error — = - ^' + -t!IL\ 

Q 2[ Ap p H J 

is thus within 2.5% for a mixture density p H = 300 kg/m 3 . 
[0095] The sampling probe 20 has an opening area 
(Aprobe), or ratio of its area to the pipe area (Ap ipe ), 
which is known accurately to within 2% and is maintained 
unchanged over a long period. To achieve this, the probe 
tip 18 is made from damage-resistant material such as 
tungsten carbide or cobalt alloy. 

[0096] Consideration is also given to the size of the 
opening of the sampling probe tip 18. In order to keep 
the relative error in the area of the opening of the 
probe, or in the area ratio to within 2%, a 1% accuracy 
is required in the diameter of the opening of sampling 
probe 20. Therefore for a small opening diameter 
sampling probe 20 (e.g. 1.5 mm), the tolerance in 
diameter has to be very high (< 15 urn) . Increasing the 
opening diameter of the sampling probe 20, for example to 
3 mm, reduces the demand on the machining tolerance for 
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the probe (to < 30 urn for the 3 mm example above) . For 
high GLR flows, a larger probe size may also be 
beneficial as it will increase the amount of liguid 
sampled and therefore ease the sizing of the Coriolis 
mass flow meters 30a and 30b for use at higher ranges of 
gas and liguid flow rates in the mini-separator 50. 
[0097] It is preferable to have a uniformly mixed wet- 
gas flow at the sampling cross section. This is achieved 
by placing a flow conditioner 200 upstream of sampling 
probe 20. Preferred embodiments of the flow conditioner 
200 are described below. 

[0098] It is also preferable that the stream at the 
Venturi portion of the Vx flow meter 40 is also 
homogenous wet-gas flow (i.e. with negligible gas-liguid 
slip), and so a further flow conditioner may be provided 
upstream of the Vx flow meter 40. 

[0099] With the on/off valve 35 open and operating in 
conditions of well-mixed wet-gas flows in which there is 
low liguid mass (i.e. Q L /Q G < 1) the total mass flow rate 
measurement Q from the Vx flow meter 40 

( Q = KyjAp.p H =Q G + Q L =Q G {l + Q L /Q G ) ) is first used to set an 
approximate 'dry' gas sample extraction mass flow rate 
Q SG = [A probe I A pjpe )x Q by adjusting flow control valve 32 in 
the mini-separator 50 through use of the actuator 34. 
[00100] The mass flow rate of the 'dry' gas sample flow 
Q SG , after gas-droplet separation in the gas-drop 
separator 51, can be directly measured using a Coriolis 
flow meter 30a. 

[00101] The sampled liguid mass flow rate Q SL is 
measured using an instrumented liguid collection system 
52 in the mini-separator 50. Precise achievement of 
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isokinetic sampling has less effect on the sample liquid 
mass flow rate Q s ,l than on the sample gas flow rate 
because of the higher inertia of liquid droplets compared 
to the gas. 

[00102] In the mini-separator 50 the sampled liquid 
mass flow rate Q SfL can be derived from the change in the 
liquid level over time h(t) in the liquid collection 
vessel 52 from a Dp level sensor (not shown) 
( p = p s L gh; Q S L = (A L / g).{dp / dt) ) , where A L is the cross-sectional 
area of the liquid collection vessel 52). 
[00103] However, in the embodiment illustrated in 
Figure 5, Q S , L is also measured using a second Coriolis 
mass flow meter 30b, which also gives the liquid mixture 
density Ps,i- For example the second Coriolis mass flow 
meter 30b may be a Promass 80A™ Coriolis meter as made by 
Endress-Hauser . This meter can measure liquid mass 
flowrate to ±0.1%, gas mass flowrate to ±0.5% and density 
to ±1 kg/m 3 , and is reported to measure flow rates as 
small as 0.1 kg/h. 

[00104] An extra linear array of MEMS (Micro Electro 
Mechanical Systems) densitometer can be installed along 
the height of the liquid collection vessel 52, allowing 
oil and water densities (p s ,o, Ps r w) to be determined from 
the separated oil-water interface profile at separator 
temperature and pressure (T s , p s ) • 

[00105] It is preferable that the liquid is separated 
at line pressure not following a pressure reduction to 
minimize mass transfer between the phases. It is also 
preferable that the liquid sample mass flow rate 
measurement system, whether it is a Coriolis meter or a 
measurement bottle has an error of 2% or less to achieve 
good results in the closed-loop feedback system described 
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below (a better than 2% accuracy is necessary to obtain 
an accurate measurement of the sampled liquid flow rate 
Q SrLl and of the main-pipe liquid rate Q L = {A pipe I A prob Jx Q s L ) . 
[00106] The water cut ( y w ) can then be determined from 
the liquid mixture density relation p SL = y w p sw + {^~7w)Ps,o • 
Alternatively the water cut can be measured by other 
means, for example using microwave or optical sensors, 
such as the RedEye™ infrared water-cut meter from 
eProduction Solutions, Texas US, prior to complete 
separation of the liquids. 

[00107] The first measure of the liquid sample mass 
flow rate Q SL is then combined with the measured sample 
gas mass flow rate in the controller 60 to obtain the 
total sample mass flow rate Q s . The controller 60 then 
compares this with the current main pipe total mass flow 
rate as measured by the main pipe flow meter 40. The 
sample extraction rate 19 is then varied as necessary by 
adjusting the flow control valve 32 in the mini-separator 
50, through use of the actuator 34, so as to cause the 
sample mass flow rate to approach [A probe I A pipe )x Q , and thus 
isokinetic sampling conditions. 

[00108] The above process can then be repeated for 
successive measurements of the total sample mass flow 
rate, thereby achieving convergence with the isokinetic 
sampling conditions. The process may be run effectively 
continuously throughout the sampling process, and thereby 
adjust the sampling flow rate as necessary when the flow 
of the main stream varies for any reason. In this 
respect, the term "continuous" includes repeated discrete 
applications of the measuring and correcting steps above 
at regular intervals, these intervals being shorter than 
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the time scale over which significant changes in the 
total mass flow rate of the main stream occur. A typical 
interval is less than 30 sec, and preferably less than 10 
or 5 sec. 

[00109] For the maintenance of isokinetic sampling in 
this system, the existence of a pressure difference 
between the main pipe 220 and the mini-separator 50 does 
not affect the sampling. Mass transfer between the 
liguid and gas phases can occur without prejudicing the 
sampling control process since the process is based on 
total mass flows and not on the volume rate of a single 
phase . 

[00110] However, it is preferable that the pressures in 
the mini-separator 50 and in the main pipe 220 are 
approximately equal, so that no mass transfers between 
the gas and liquid phases occur between the pipe and the 
separator. In this situation, the main pipe gas and 
liquid mass flow rates can be obtained as follows: 
[00111] Total mass flow rate Q = K^Ap.p H , obtained from 
the Vx Venturi Ap measurement and the high-energy mixture 
density p H measurement; 

[00112] Liquid mass flow rate Q L ={A pipe I A probe )xQ SL , 
obtained from the isokinetic sampling; 

[00113] Gas mass flow rate Q G =Q-Q L =Q- [A pipe I A probe )x Q SL , 
preferably obtained from the Vx total flow measurement 
and isokinetic sampling, also obtainable purely from the 
isokinetic sampling as Q G = (A pipe I A prohe )x Q s G , which can 
therefore be used as a quality control check. 
[00114] As a result, the PVT relations are no longer 
reguired to convert the liquid mass flow rate measured at 
separator conditions (p s , T s ) to that at main stream 
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conditions (p, T) in order to obtain gas and liquid mass 
flow rates from the Vx measurements and interpretation 
models. However, PVT relations are required to convert 
the mass flow rates at line conditions (p, T) to volume 
flow rates at standard conditions. 

[00115] The mini-separator 50 also provides the 
following measurements of the sampled flow: 
[00116] Densities p s , G , Ps,l, Ps,o, Ps,w (of gas, liquid, 
oil and water samples respectively) obtained from the 
Coriolis gas and liquid flow meters 30a and 30b and the 
MEMS densitometer array reading from the liquid 
collection vessel 52; 

[00117] Water cut 7 W ={p s ,L~ Ps,o)^iPs,w ~ Ps,o) calculated from 
the densities above; 

[00118] GLR = (Q sg I ' p StG )/{Q StL / 'p StL ) calculated from the 

separator gas and liquid volume flow rates; 

[00119] Liquid mass loading m = Q SL IQ SG calculated from 

the separator flow rates - ideally m < 1 for a well mixed 

wet-gas regime which the system is being used to analyse; 

and 

[00120] PVT data from the samples obtained from the 
separated liquid flow and separated (i.e. 'dry') gas flow 
in the mini-separator 50. 

[00121] Another embodiment of the apparatus and method 
of the present invention is illustrated schematically in 
and with reference to Figure 6 . 

[00122] As with the embodiments described above, it is 
preferable that the wet-gas flow at the sampling cross- 
section is well-mixed, and this is achieved using mixer 
or flow conditioner 200 positioned upstream of the 



57.0503 U.S. PCT 



26 

sampling cross-section. Again, preferred embodiments of 
the flow conditioner are described below. 
[00123] An averaging Pitot-tube (APT) 44 is used to 
derive the main pipe total mass flow rate when operating 
at impact-pressure measuring (Pitot) mode, as well as as 
a sampling probe. The APT 44 has multiple holes, for 
example four or six holes, specifically arranged along 
the diametrically aligned tubing. For example, for a 
six-hole probe the holes are preferably placed on the 
probe tubing at the following positions: -.913R, -.707R, 
-.408R, +.408R, +.707R, +.913R, where R is the radius of 
the pipe. The hole positions preferably correspond to 
the centres of the sampled egual area annuli of the pipe 
circular cross section. As with the previous embodiment, 
for sampling accuracy it is preferable that the area of 
the hole openings of the APT 44 (A probe = N hole x A hole ) , 
relative to the pipe area A pipe at the sampling cross- 
section, is known to within 2% and does not change over 
the time of testing and sampling. This can be achieved 
by making the APT 44 from damage resistant material such 
as tungsten carbide or cobalt alloy. 

[00124] One example of an APT is the Rosemount Probar™ 
insertion meter operating at impact-pressure measuring 
(Pitot) mode only, which claims a relative accuracy of 
±1.1% in volumetric flow rate measurements for steam, 
liguid and gas . APTs can also be obtained from 
Thermocouple Instruments Ltd. Of Cardiff, UK. 
[00125] An APT probe has very low pressure loss and 
generally costs less than a Venturi flow meter. The pipe 
averaging characteristics of the APT 44 reduce any 
variations in the cross-pipe flow at the sampling cross 
section . 
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[00126] A high energy gamma-ray densitometer 43 is 
installed upstream of the APT probe 44 and downstream of 
the flow conditioner 200 to obtain the density p H . The 
densitometer 43 can have an accuracy within 10 kg/m 3 from 
a 30 mCi 137 Cs source (with high energy gamma-ray at 
661.66 KeV) and a Nal scintillation detector. 
[00127] The APT 44 is first set to operate in normal 
impact pressure measurement mode, with on/off valve 35 
set in the off position and with no sampling flow being 
taken. The APT measures the differential pressure Ap 
across the main-pipe static pressure tap and the APT 
probe static pressure tap (typically to within 1% 
accuracy) , and the main pipe total mass flow rate is 
calculated as Q = K^Ap.Pu (here K 2 is a calibration 
factor dependent on the probe geometry and flow Reynolds 
number) . 

[00128] The APT 44 is then switched to operate in 
sampling mode by opening the on/off valve 35 and by 
adjusting flow control valve 32 through actuator 34. The 
averaging nature of the APT 44 means that the main flow 
stream is sampled at multiple points across the diameter 
of the sampling cross section. Therefore variations in 
the cross-pipe GLR profile, which may be found downstream 
of an imperfect flow conditioner, can be largely 
overcome. Alternatively, when an APT is used, as in the 
current embodiment, the performance reguirements of the 
flow conditioner can be made less stringent. 
[00129] Operating in conditions of wet-gas flows in 
which there is low liguid mass (i.e. Q l /Qg < 1) the total 
mass flow rate measurement Q from the gamma-ray 
densitometer 43 and the APT 44 operating in impact 
pressure measurement mode ( 
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Q = KjAp.p H =Q G +Q L =Q G {l + Q L /Q a ) ) is first used to set an 
approximate 'dry' gas sample extraction mass flow rate 
Qs,g 0 = {A pro be / A pipe )xQ by adjusting flow control valve 32 in 
the mini-separator 50. 

[00130] The remainder of the sampling process is the 
same as that described in relation to Figure 5 above. In 
particular the flow control valve 32 is adjusted via the 
actuator 34 so that the sampling flow 19 is substantially 
isokinetic with the main pipe flow stream 11. 
[00131] In a further embodiment, a single-hole sampling 
probe 20, such as that used in the embodiment described 
in relation to Figure 5 above, can also be set to operate 
in impact pressure mode by switching valve 35 off and 
measuring the differential pressure Ap across a main pipe 
static pressure tap and a static pressure tap on the 
probe 20. This measurement of Ap could then be used to 
derive an independent estimate of the main pipe total 
flow rate using the density p H measured by densitometer 
42 . 

[00132] In a further embodiment, in applications where 
the probe pressure loss (i.e. from the probe opening to 
the probe pressure tap) is relatively small, the 
isokinetic sampling condition can be substantially 
achieved by nulling the Ap described above by adjusting 
the sampled flow rate of the sampled portion through the 
use of flow control valve 32. In this embodiment, the Ap 
signal can be used as a control signal to control 
actuator 34, thus avoiding the use of flow meters 30a and 
30b. Note that this embodiment can be adapted to either 
the single hole probe shown in Figure 5 or the multiple 
hole pitot tube probe shown in Figure 6 . 
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[00133] Note that the multiple-hole pitot tube sampling 
embodiments described in connection with Figure 6 can 
also be used in applications where the independent 
measurement is made using Vx type meters as described in 
connection with Figure 5 . 

[00134] Figure 7 is a diagram of a pipe section which 
is arranged to provide a relatively uniform cross-pipe 
profile at the sampling cross-section, particularly with 
regard to the GLR. The pipe section shown has a flow 
conditioner 200 for providing such a cross-section. 
[00135] The flow conditioner 200 has three elements: a 
swirl inducing portion, which is provided in this 
embodiment by the tangential flow input 201; a gas flow 
core straightener 202; and an orifice plate 203. 
[00136] The tangential flow input 201 causes the entire 
flow to swirl, thereby creating a substantially uniformly 
circumf erentially distributed liquid film on the pipe 
wall 205. This removes variations due to the inflow 
conditions, such as non-axisymmetr ic gas and liquid flow 
distributions . 

[00137] The core of the flow, which is mostly gas, is 
then straightened by the core straightener 202 positioned 
downstream of the tangential input 201. An example of 
such a core straightener 202 is illustrated in Figure 9 
and comprises a conventional straightener 212 (such as 
those described in US 3840051, US 5309946 and US 5762107) 
positioned with an annular gap 213 between the 
straightener 212 and the pipe wall 205 so as not to 
interfere with the liquid flow along the pipe wall. The 
core straightener 202 removes the swirl in the central 
gas core of the flow, thereby reducing the amount of 
liguid re-deposited on the wall after the orifice plate 
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mixing. This facilitates the obtaining of consistent 
gas-liguid flux profiles at the sampling cross-section. 
[00138] Downstream of the core straightener 202, an 
orifice plate 203 is positioned, which causes the liguid 
film flowing on the pipe wall to be shed back into the 
now straightened gas core. The diameter (5-ratio (the 
ratio of the open diameter of the orifice plate 203 to 
the diameter of the pipe) is preferably between 0.85 and 
0.95. A further orifice plate (not shown) may be 
positioned downstream of orifice plate 203, and sheds any 
liguid film re-deposited on the wall after the orifice 
plate 203. 

[00139] The liguid shed from the pipe wall by the 
orifice plate 203 is then subjected to an acceleration- 
induced catastrophic break-up by the gas flow which has a 
high axial velocity. In the region immediately 
downstream of the orifice plate 203, a liquid rich cone 
is formed at the centre of the flow. Downstream of this 
cone, the turbulent kinetic energy is high and hence the 
secondary turbulent break-up controls the droplet size 
distribution. Further downstream, as the turbulent 
kinetic energy decays quickly, turbulent break-up ends 
and coalescence effects determine the liquid droplet 
size. Thus the cross-pipe liquid flux distribution 
downstream of the orifice plate 203 is totally different 
from that upstream of the plate 203. 

[00140] The sampling probe is preferably placed at an 
axial distance L downstream of the orifice plate 203, the 
distance L being chosen so as to avoid sampling at the 
liguid-rich core region caused by the shedding of the 
liguid film by the orifice plate 203, and to allow for 
some break-up due to turbulence and diffusion, leading to 
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smaller drop sizes and more uniform cross-pipe GLR. The 
preferred values of L are ID to 3D. 

[00141] The use of an orifice-plate type mixer 203 is 
preferable to the flow-focusing type mixers, such as 
those disclosed in GB 2301297 and US 5894080. An easy- 
to-machine orifice-plate mixer 203 serves to shed the 
liquid film from the pipe wall, which is then subjected 
to acceleration-induced catastrophic break-up by the 
incoming high axial-velocity gas flow. For a 4" pipe 
diameter (D = 100 mm) down-flow example at low pressures, 
a much more uniform cross-pipe GLR has been 
experimentally measured 3D downstream of a single 
orifice-plate 203 with diameter ratio (3= 0.95 (see Figure 
2b) than that obtained with the flow-focusing mixer 15 
shown in Figure 2a (which is similar to that described in 
GB 2301297) . As a result, the GLR error at the pipe 
centre is reduced by about five-fold over the flow- 
focusing type mixer shown in Figure 2a. The orifice 
plate and the flow-focusing mixer were both tested at the 
same CGR of 60 bbl/MMscf, at a gas rate of 1.2 MMscf/day. 
[00142] The embodiment of Figure 7 shows the flow 
conditioner and sampling probe 20 implemented in a pipe 
section 210 of smaller diameter than the main pipe 220 
(for example a section of 2 inch (D = 50 mm) rather than 
4 inch (100 mm) diameter) . For a given mass flow rate of 
a well stream, deploying the sampling probe 20 at a small 
diameter pipe section 210 achieves the highest gas flow 
momentum (p G U G 2 ) and hence the highest gas Weber number 
We D = p G U G 2 D/o (where o is the surface tension) . This may 
be beneficial as it: increases the shear for the liguid- 
film primary break-up from the wall, hence increasing 
drop formation; causes generation of small size droplets, 
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thereby reducing the drop Stokes number and gas-drop 
velocity slip; increases the cross-pipe uniformity and 
the percentage of liguid entrainment. 
[00143] An alternative embodiment of the flow 
conditioner 200 is shown in Figure 8. This embodiment is 
the same as that in Figure 6, except that the tangential 
input 201 has been replaced by a bend inlet 221 followed 
by a helical insert 211 to provide the initial swirl. 
[00144] Alternative embodiments which cause 
axisymmetric distribution of the liguid film on the pipe 
wall, by causing swirl in other ways may also be used. 
[00145] The above embodiments may form part of an 
overall sampling apparatus, being combined with the 
apparatuses and methods of earlier embodiments to provide 
a uniform cross-pipe GLR for sampling. 
[00146] While the invention has been described in 
conjunction with the exemplary embodiments described 
above, many eguivalent modifications and variations will 
be apparent to those skilled in the art when given this 
disclosure. Accordingly, the exemplary embodiments of 
the invention set forth above are considered to be 
illustrative and not limiting. Various changes to the 
described embodiments may be made without departing from 
the spirit and scope of the invention. 



